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PREFACE 

This  Memorandum  is  part  of  RAND's  continuing  VELA  Analysis 
study  for  the  Advanced  Research  Projects  Agency  to  assess  the  rel¬ 
ative  capabilities  of  various  methods  of  detection  of  nuclear  explo¬ 
sions  above  the  earth's  surface.  The  riometer  technique  is  one 
possible  method  for  the  detection  of  nuclear  explosions  at  high 
altitude  and  in  deep  space.  The  maximum  range  or  altitude  at  which 
a  1-kT  X-ray  impulse  can  be  detected  by  measuring  changes  in  high- 
frequency  cosmic  noise  absorption  utilizing  a  riometer  is  determined 
in  the  Memorandum.  It  is  shown  that  the  detection  range  is  a  function 
of  the  riometer  time  constant. 

This  Memorandum  was  submitted  to  Dr.  Robert  Frosch,  Deputy 
Director  of  the  Advanced  Research  Projects  Agency,  for  inclusion  in 
the  IEEE  special  issue  on  nuclear  detection. 


SUMMARY 


This  Memorandum  considers  the  detection  of  overhead  nuclear 

explosions  in  space  by  riometer  measurement  of  the  effects  of  explosion- 

induced  ionization  on  the  absorption  of  high-frequency  cosmic  noise. 

The  induced  ionization  is  calculated  as  a  function  of  height  and  time 

for  a  1-kT  X-ray  impulse  having  temperatures  of  .5,  1.0,  and  2.0  kev 

35 

and  having  burst  heights  from  10  to  10  km.  The  absorption  at  each 

height  due  to  the  induced  ionization  is  then  calculated  at  a  radio 

frequency  of  30  Me.  In  order  to  determine  the  change  in  the  effective 
antenna  temperature  as  a  function  of  time  at  the  riometer,  the  total 

absorption  through  the  ionosphere  is  then  determined.  The  response 

of  a  riometer  with  several  time  constants  is  found  for  this  sharp 
change  in  received  cosmic  noise  power.  Detection  ranges  are  calculated 
for  each  X-ray  temperature  and  several  possible  riometer  time  constants 
based  on  the  criterion  that  the  nuclear  burst  must  produce  a  1-db 
peak  response  in  the  diode  current  of  the  riometer.  These  detection 
ranges  are  tabulated  and  compared  with  those  that  were  found  neglecting 
riometer  response  time. 
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I.  INTRODUCTION 

The  riometer^  is  a  very  stable  receiver  used  to  measure  the 
relative  change  in  cosmic  noise  arising  primarily  from  synchrotron 
radiation  within  our  own  galaxy.  The  relative  changes  in  received 
noise  power  are  measured  quite  accurately  by  continuously  comparing 
the  antenna  noise  temperature  with  tne  noise  temperature  of  a  controlled 
noise  source.  In  the  absence  of  the  earth's  ionosphere,  the  intensity 
of  the  cosmic  radio  noise  from  any  given  large  sector  of  the  galaxy 
is  constant  with  only  the  usual  statistical  fluctuations.  In  the 
presence  of  the  ionosphere,  the  portion  of  the  noise  spectrum  below 
f^F2  is  basically  cut  off,  while  the  frequency  components  higher  than 
fQF2  pass  through  the  ionosphere  and  are  attenuated.  The  amount  of 
attenuation  introduced  by  the  ionosphere  is  frequency  dependent  and 
varies  in  time  with  the  amount  of  ionization  in  the  ionosphere.  When 
a  high-altitude  nuclear  burst  occurs,  the  ambient  ionization  in  the 
ionosphere  may  be  enhanced,  causing  increased  attenuation  of  the  cosmic 
radiation  that  may  be  detected  by  a  riometer  on  the  earth. 

A  large  body  of  literature  is  available  which  calculates  and 

estimates  the  effects  of  explosion-produced  ionization  on  radio  waves.  A 

few  of  these  papers  which  have  been  used  extensively  for  background 

(2) 

in  this  paper  are  mentioned  below.  Latter  and  LeLevier  have  cal¬ 
culated  the  vertical  absorption  applicable  to  the  study  of  riometer 
measurements  due  to  X-ray,  y-ray,  and  p-ray  radiation  of  high-altitude 

(3) 

nuclear  bursts  neglecting  any  ambient  ionization.  Crain  and  Booker 
developed  a  simple,  but  general,  relation  for  the  prompt  vertical 
absorption  illustrating  its  dependence  on  operating  frequency,  angle 
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of  incidence,  maximum  induced  ionization,  and  the  X-ray  radiation 

temperature.  The  time  dependence  of  the  ionization  impulse  decay 

(4) 

was  also  discussed.  This  work  was  extended  by  Crain'  in  a  paper 
which  pointed  out  the  nonlinear  effects  in  the  rate  of  ioniza¬ 
tion  decay  that  may  also  occur  depending  on  the  intensity  of  the 
X-ray  impulse.  An  analytical  expression  for  the  decay  of  the 
ionization  for  the  special  case  of  equal  electron- ion  and  ion- ion 
recombination  rates  was  developed  by  LeLevier.^ 

In  this  Memorandum,  the  X-ray  induced  ionization  is  calculated  for  a 
1-kT  X-ray  yield  for  a  range  of  altitudes  and  hypothetical  X-ray  tempera¬ 
tures.  The  decay  of  this  ionization  as  a  function  of  time  is  then 
computed,  as  well  as  the  resulting  absorption  as  a  function  of  height. 

The  total  one-way  absorption  through  the  ionosphere  is  also  calculated 
as  a  function  of  time  after  the  burst,  and  the  rlometer  response  to 
this  attenuation  is  computed  taking  into  account  its  hypothetical 
time  constant  a  id  antenna  pattern. 

From  the  calculation  of  riometer  response  to  an  X-ray  impulse', 
detection  ranges  are  calculated  which  are  used  to  draw  conclusions 
as  to  the  utility  of  the  riometer  in  detecting  nuclear  explosions  in 
space.  Explosions  at  altitudes  below  about  1000  lan  are  not  considered. 
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II.  THEORETICAL  X-RAY  IONIZATION  IMPULSES 
(23 

Standard  procedures'  7  exist  for  obtaining  the  energy  deposited 

at  each  altitude  due  to  a  given  flux  of  X-ray  radiation  from  a  nuclear 

explosion.  Using  these  procedures,  the  vertical  ionization  depositions 

3  5 

were  calculated  for  a  1-kT  X-ray  yield  at  altitudes  from  10  to  10 
km  which  have  (hypothetical  X-ray)  radiation  temperatures  of  .5,  1.0, 
and  2.0  kev. 

The  computed  values  of  the  electron  density,  AN,  at  vertical 

incidence  versus  height  for  a  1-kT  X-ray  impulse  at  a  height  of  1000 

km  at  the  instant  after  the  impulse  occurs  are  shown  in  Fig.  1  for 

three  X-ray  temperatures.  The  curves  of  Fig.  1  will  hold  almost 

equally  well  for  a  100-kT  burst  at  10,000  km  since  the  initial,  or 

(33  2 

prompt,  ionization  scales'  7  as  Y/R  where  Y  is  the  yield  in  kT  and 
R  is  the  distance  from  the  burst  altitude  to  the  altitude  of  interest 
in  km.  The  initial  ionization  curves  of  Fig.  1  hold  equally  well 
for  day  or  night. 

The  curves  of  Fig.  1  demonstrate  that  the  cooler  the  radiation, 
the  higher  the  ionization  is  deposited.  To  clarify  the  magnitudes 
of  the  ionization  impulses  under  discussion,  Table  1  is  given  to 
illustrate  the  peak  electron  density  for  each  height  and  temperature 
of  the  hypothetical  1-kT  X-ray  impulse. 
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Table  1 

APPROXIMATE  PEAK  VALUES  OF  1-kT  X-RAY  IMPULSES 


1 

Burst  height 
(km) 

AN  (elec/cc) 

Temperature  (kev) 

0.5 

(91  km)  ^a,) 

1.0 

(82  km)(a) 

2.0 

(71  km/a) 

io3 

2.9  x  106 

! 

2.8  x  106 

1 

2.3  x  IO6 

3  x  103 

2.9  x  105 

2.8  x  IO5 

2.3  x  IO5 

io4 

2.4  x  104 

2.4  x  104 

2.0  x  IO4 

u 

3  x  10 

3 

2.6  x  10 

2.6  x  IO3 

2.2  x  104 

io5 

2.4  x  IO2 

2.4  x  IO2 

2  x  IO2 

(a) 


Height  of  peak 
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III.  DECAY  OF  AN  IONIZATION  IMPULSE 

(3) 

Using  the  standard  rate  equations'  for  atmospheric  ionization 

and  commonly  used  values  of  the  recombination  coefficients^,  the 

decay  of  X-ray  deposited  ionization  was  computed  versus  time  for 

heights  from  30-140  km.  The  values  obtained  are  typical,  and  the 

final  result  is  relatively  independent  of  reasonable  changes  in  the 

(4) 

recombination  coefficients  such  as  suggested  by  Crain  *  or  recent 
laboratory  measurements . ^ ^  The  collision  frequency  as  a  function 

/g\ 

of  height  was  taken  from  a  paper  by  Nicolet'  *  (the  lower  of  his 

estimated  values  for  each  altitude  above  85-km  altitude). 

The  set  of  first-order  differential  rate  equations  was  solved 

numerically  on  the  computer  using  an  Adams-Moulton  variable  step  size 

predictor-corrector  procedure.  The  solutions  led  to  the  distribution 

of  enhanced  ionization  as  a  function  of  height  and  time  for  the 

hypothetical  nuclear  bursts  of  interest.  Figure  2  illustrates  the 

case  which  demonstrated  the  most  dramatic  change  with  time.  In 

Fig.  2  it  is  seen  that  the  electrons  deposited  at  t  *  0  are  quickly 

attached  to  the  neutral  particles  at  altitudes  below  65  km.  A  detailed 

discussion  of  the  decay  results  obtained  is  beyond  the  scope  of  this 

Memorandum.  Readers  interested  in  the  details  of  decay  and  persistance  of 

(4) 

an  impulse  of  ionization  are  referred  to  a  paper  by  Crain.  ' 


X-ray  yield  =  1  kT 

Burst  height  =  1000  km 

Radiation  temperature  =  2  kev 
Nighttime 


Fig. 2— Decay  of  ionization  deposition  versus  altitude  and  time 
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IV.  COSMIC  NOISE  ABSORPTION 

The  impulses  of  ionization  produced  by  a  high-altitude  X-ray 
burst  will  cause  increased  absorption  at  a  given  altitude  according 
to  the  absorption  expression 


4.6  x  10 


(N  +  AN)v 
a 


■)  -  —  •)  1 

'  V  ♦  -,2'  1 


where 

A'  m  added  absorption  in  db  per  km 

*  ambient  electron  density  (elec/cm') 
v  ■  collision  frequency  (coll/sec) 

0)  ■  angular  operating  frequency  (rad/sec) 

■  real  part  of  the  ambient  refractive  index 

■  value  of  ii  after  the  impulse 
AN  ■  X-ray  deposited  ionization 

Figures  3  and  4  illustrate  the  absorption  value  A1  as  a  function 

4 

of  height  and  time  for  a  1-kT,  .5  kev  X-ray  impulse  at  10  km  for 
day  and  night  respectively. 

The  one-way  absorption,  A,  of  galactic  noise  is  then  given  simply 

by  140  km 

A  -  [  A'  dh  (2) 

J  30  km 


Values  of  A  as  a  function  of  time,  A(t) ,  were  calculated  for  a  1-kT 

3  3  4  4 

X-ray  impulse  occurring  at  altitudes  of  10  ,  3  x  10  ,  10  ,  3  x  10  , 

and  10^  km.  They  were  calculated  for  temperatures  of  .5,  1.0,  and 

2.0  kev  for  both  day  and  night,  making  a  total  of  30  cases  studied. 

The  values  of  A  were  examined  primarily  for  an  operating  frequency,  f  ■  30 


40  60  80  100  120  140 

Height  (km) 

Fig.3 — Daytime  absorption  at  30  Me  versus  height  and  time 
for  a  decaying  ionization  deposition 


(db/km ) 


40  60  80  100  120 


Height  (km) 

Fig. 4— Nighttime  absorption  at  30  Me  versus  height  and  time 
for  a  decaying  ionization  deposition 
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Examples  of  the  results  obtained  for  A(t)  are  given  in  Fig.  5  for 
daytime  conditions  and  in  Fig.  6  for  the  nighttime.  Figure  5  shows 
how  the  daytime  one-way  vertical  absorption  quickly  decays  to  an 
almost  constant  value  in  the  interval  of  interest  due  to  the  existence 
of  photoattachment;  while  the  nighttime  decay,  in  the  absence  of  photo¬ 
attachment,  has  a  more  negative  slope  of  dA/dt. 

An  Interesting  feature  of  Fig.  6  is  that  after  approximately  2.2 
sec,  the  value  of  A  is  higher  for  the  1.0  kev  X-ray  impulse  tempera¬ 
ture  than  for  the  2.0  kev  impulse.  This  is  due  to  the  fact  that  the 
higher  temperature  deposits  more  ionization  at  lower  altitudes  where 
the  attachment  rate  is  high.  This  does  not  occur  during  the  day  since 
the  photodetachment  causes  slightly  more  presistance  at  the  lower 
altitudes  for  times  which  are  short  compared  to  the  recombination 
relaxation  time. 

The  nonlinear  recombination  decay  effect  which  occurs  when 

(4) 

initial  depositions  are  quite  intense  has  been  demonstrated  by  Crain. 
Figure  7  illustrates  this  effect.  Note  how  the  low-altitude  impulse 
decays  at  a  higher  rate  than  the  higher-altitude  impulse  due  to  the 
higher  initial  ionization,  as  shown  in  Table  1. 


8 


0 


2  4  6 


Time  after  impulse  (sec) 


Fig. 5— Daytime  one-way  vertical  absorption  at  30  Me  versus 
time  due  to  the  decay  of  X-ray  induced  ionization 


Bunt  height  =  10,000  km 
Yield  =  1  kT 


1  key 

2  kev 


0.5  kev 


Fig. 6— Nighttime  one-way  vertical  absorption  at  30  Me  versus 
time  due  to  the  decay  of  X-ray  induced  ionization 
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t  (sec) 

Fig. 7— One-way  vertical  absorption  at  30  Me  versus  time 
for  a  1  kT,  0.5  kev  nighttime  X-ray  impulse  at  several  altitudes 
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V.  RIOMETER  RESPONSE  TO  ABSORPTION  IMPULSES 


A  typical 


rioraeter  operates  by  holding  the  equality 


(9) 


T  =  T  +  5800  IR 
A  R 


(3) 


where 

=  riometer  antenna  temperature,  °K 

T  =  servo  diode  resistor  temperature 
K 

I  =  servo  diode  current 

R  =  value  in  ohms  of  the  servo  diode  resistor 
The  value  5800  involves  the  electronic  charge,  e,  and  Boltzmann's 
constant,  Relative  changes  in  antenna  noise  power,  AP  =  k  T  B 

A  AW 

(where  By  is  the  bandwidth) ,  cause  proportional  changes  in  the  servo 
diode  current  so  that 


AP 

10  logio  (~p^)  “  10  logio  (^i)  ;  steady  state 


(4) 


Th'’  above  equality,  however,  does  not  hold  instantanteously  since 

riometers  require  integration  over  finite  time  in  order  to  suppress 

the  statistical  fluctuations  in  P  .  The  upper  bound  of  the  variance 

A 

of  the  riometer  response  is^^ 


where 


0^  =  variance  of  the  diode  current 
I  =  average  diode  current 
By  *  predetection  bandwidth 
Tc  =  postdetection  integration  time  constant 


(5) 
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If  a  predetection  bandwidth  of  10,000  cycles  per  sec  is  used,  it  is 
found  from  the  above  relationship  that  the  integration  time  constant 
must  be  1  sec  for  the  variance  to  remain  at  a  1  per  cent  level.  This 
is  why  most  integration  time  constants  used  in  riometers  range  from 
.2  to  27  sec.  The  bandwidth  cannot  be  greatly  increased  due  to 
interference  problems;  therefore,  the  riometer  must  have  integration 
time  constants  that  are  rather  long  compared  to  the  time  constants  of 
the  receiver  AGC  circuit. 

The  transfer  function  of  a  riometer  must  be  of  the  nature 

H(t)  =  (^-  e‘t/Tc)  (6) 

c 

where 

H(t)  =  transfer  function 

T  =  integration  time  constant  (sec) 
c 

t  =  time  (sec) 

The  response  of  the  riometer  is  then 

AI(t)  *  K  J  AP(t)  e'(t  '  t)/Tc  dr  (7) 

c  J  o 

where 

AI(t)  =  change  in  the  servo  diode  current 
K  =  a  proportionality  constant 
AP(t)  =  change  in  noise  power  received 
T  =  a  dummy  variable  in  time 

t  =  time  measured  from  the  instant  the  change  in  AP  occurs 
In  the  calculation  of  change  in  AI(t),  the  proportionality  constant 
was  unimportant  since  only  relative  changes  of  —■  (in  db)  were  of  interest. 


Figure  8  presents  an  example  of  the  rlometer  response  to 
X-ray  Impulse  Induced  absorption  phenomena  for  two  different  integra¬ 
tion  time  constants.  In  Fig.  8,  it  is  seen  that  the  shorter  the 
integration  time  constant,  the  higher  the  peak  response  of  the 
rlometer.  However,  as  discussed  previously,  the  time  constant  can¬ 
not  be  shortened  without  bound  due  to  the  statistical  nature  of  the 
cosmic  radio  noise  and  the  increased  possibility  of  false  alarms. 
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(qp)  a6uoip  4»N 


t  (sec) 

Fig. 8— Typical  riometer  response  to  an  X-ray  absorption  event  versus 
time  for  two  integration  time  constants 
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VI.  MAXIMUM  DETECTION  RANGE  OF  X-RAY  IMPULSES 

(2) 

Latter  and  LeLevier  reported  values  of  R  (where  R  is  the  dis¬ 
tance  in  km)  at  which  a  1-kT  yield  of  X-rays  produces  a  one-way 
vertical  absorption  of  1  db  at  an  operating  frequency  of  30  Me. 

These  values  of  R  were  reported  for  X-ray  temperatures  of  2,  1,  .5, 
and  .1  kev  as  well  as  for  several  angles  of  incidence  of  the  X-ray 
energy  to  the  vertical.  These  values  of  R  were  based  strictly  upon 
the  initial  or  prompt  absorption. 

In  utilizing  riometers  to  measure  the  cosmic  noise  absorption, 
the  initial  absorption  cannot  be  dealt  with  since  the  riometer  will 
not  respond  instantaneously  to  the  decreased  power  at  its  antenna 
terminals.  Therefore,  utilizing  the  previously  discussed  work  of  the 
absorption  versus  time  and  the  convolved  riometer  response  to  this 
absorption,  a  new  set  of  values  for  R  is  reported  in  Table  2  for  vertica 
incidence  only  but  as  a  function  of  several  possible  time  constants. 

The  results  were  obtained  using  peak  value  of  plots  similar  to  those 
shown  in  Fig.  8  for  both  daytime  and  nighttime  conditions. 

Table  2 

VALUES  OF  R  AT  30  Me  FOR  WHICH  A  1-kT  YIELD 
_ CAUSES  A  1-db  PEAK  RIOMETER  RESPONSE 
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It  is  seen  from  Table  2  that  the  maximum  distance  at  which  a 

1-kT  X-ray  impulse  can  be  detected  in  this  temperature  range  using 

cosmic  noise  absorption  is  slightly  over  3  earth  radii.  The  daytime 

and  nighttime  detection  ranges  are  comparable,  although  the  range  is 

always  slightly  less  at  night. 

The  small  differences  in  Table  2  between  the  T  =  0  values 

c 

deduced  in  this  paper  and  those  given  in  Ref.  2  arise  from  the 
use  of  somewhat  different  collision  frequency  profiles  and  methods 
of  integration. 
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VII.  PRACTICAL  CONSIDERATIONS  IN  THE  UTILIZATION  OF  RIOMETERS 

Many  practical  considerations  exist  in  the  utilization  of  rio- 

meters  as  part  of  a  system  for  detecting  nuclear  explosions  in  addition 

to  their  limited  detection  capability.  Some  of  these  are:  limited 

visibility,  since  most  absorption  occurs  in  the  D  region  at  30  Me; 

other  sources  of  noise,  such  as  man-made  and  lightning  pulses;  and 

natural  absorption  events  due  to  solar  flares  which  may  cause  false 

alarms.  Extensive  studies  by  Dyce  et  al.^^  have  been  carried  out 

to  examine  these  practical  considerations. 

Utilization  of  the  proper  antenna  pattern  with  low  side  lobes 

is  quite  important  in  order  to  average  over  a  large  part  of  the  sky 

and  yet  discriminate  against  man-made  interference  arriving  at  low 

elevation  angles.  The  pattern  of  the  antenna  can  actually  change 

the  values  of  R  obtained  in  Fig.  8  since  oblique  angles  are  utilized 

in  averaging  over  the  antenna  pattern.  This  nonlinear  correction 

(12)  * 

has  been  studied  by  Weir'  '  and  Bedrosian  and  would  increase  R  by 

2 

a  factor  of  the  order  of  / 1 .5  for  a  cos  0  (0  *  aspect  angle)  antenna 

pattern. 

* 

Private  communication  from  E.  Bedrosian,  The  RAND  Corporation. 
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VIII .  CONCLUDING  REMARKS 

The  riometer  is  a  relatively  inexpensive  device  that  can  be 

( 13) 

distributed  over  the  entire  earth  for  the  measurement  of  natural 

ionospheric  absorption.  These  instruments  can  then  also  be  used  to 

detect  nuclear  bursts  in  space  as  well  as  to  provide  one  input  to 

the  picture  of  world-wide  ionospheric  conditions. 

As  a  device  for  detecting  deep  space  X-ray  impulses,  however, 

the  riometer  is  rather  insensitive  when  compared  to  other  techniques 

such  as  phase-path  sounders  or  HF-forward  propagation,  as  illustrated 

(3) 

by  Crain  and  Booker. 

While  only  the  detection  of  X-ray  impulses  in  space  is  considered 

in  this  Memorandum,  the  riometer  has  further  utility  in  detecting 

ionization  from  other  sources  associated  with  nuclear  explosions, 

(2) 

such  as  discussed  by  Latter  and  LeLevier. 
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IO.  ABSTRACT 

A  consideration  of  the  detection  of  over¬ 
head  nuclear  explosions  in  space  by  rioneter 
measurement  of  the  effects  of  explosion- 
induced  ionization  on  the  absorption  of 
high-frequency  cosmic  noise.  The  induced 
ionization  is  calculated  as  a  function  of 
height  and  time  for  a  1-kt  X-ray  impulse 
having  temperatures  of  .5,  1.0,  and  2.0 
kev  and  having  burst  heights  from  10 3  to 
10^  km.  The  absorption  at  each  height  due 
to  the  induced  ionization  is  then  calcu¬ 
lated  at  a  radio  frequency  of  30  Me.  In 
order  to  determine  the  change  in  the  effec¬ 
tive  antenna  temperature  as  a  function  of 
time  at  the  rioneter,  the  total  absorption 
through  the  ionosphere  is  then  determined. 
The  response  of  a  riometer  with  several  tine 
constants  is  found  for  this  sharp  change  in 
received  cosmic  noise  power.  Detection 
ranges  are  calculated  for  each  X-ray  tem¬ 
perature  and  several  possible  riometer  time 
constants,  based  on  the  criterion  that  the 
nuclear  burst  must  produce  a  1-db  peak  re¬ 
sponse  in  the  diode  current  of  the  rioneter. 
These  detection  ranges  are  tabulated  and 
compared  with  those  that  were  found  neglect¬ 
ing  riometer  response  tine. 
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